vectpr
vectpr(
q
m,
q
l
)
Func
Local a,
q
v
ker(
q
m-
q
l
[1]*identity(rowDim(
q
m)))
à
q
v
If dim(
q
l
)=1:Goto z
For a,2,dim(
q
l
)
augment(
q
v,ker(
q
m-
q
l
[a]*identity(rowDim(
q
m))))
à
q
v
EndFor
Lbl z
q
v
EndFunc
